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Résumé - Nous présentons dans ce papier deux nouveaux modèles de texture orientée, basés sur une nouvelle classe de champs gaussiens, 
appelés champs browniens fractionnaires localement anisotropes, à orientation locale prescrite en chaque point. Ces champs aléatoires sont une 
version locale d’une classe spécifique de champs anisotropes autosimilaires à incréments stationnaires. La simulation de telles textures s’appuie 
sur un nouvel algorithme couplant l’utilisation de la notion de champs tangents à la méthode de Cholesky ou plus récemment la méthode 
des bandes tournantes, cette dernière ayant prouvé son efficacité à générer des textures stationnaires anisotropes. Des applications numériques 
illustrent la faculté de la méthode à synthétiser des textures à orientation locale prescrite. 

Abstract - This paper présents two new models of oriented texture, based on a new class of Gaussian fields, called locally anisotropic fractional 
Brownian fields, with prescribed local orientation at any point. These fields are a local version of a spécifie class of anisotropic self-similar 
Gaussian fields with stationary incréments. The simulation of such textures is obtained using a new algorithm mixing the tangent field formulation 
with the Cholesky method or the turning band method, this latter method having proved its efficiency for generating stationary anisotropic 
textures. Numerical experiments show the ability of the method for synthesis of textures with prescribed local orientation. 


1 Textures et champs gaussiens autosi¬ 
milaires anisotropes 

La modélisation de texture est un problème difficile du traite¬ 
ment d’image. Il existe une variété de méthodes associées à dif¬ 
férents champs des mathématiques : des méthodes dites struc¬ 
turelles, statistiques, ou encore à base de patchs, etc. Nous nous 
intéressons ici aux textures aléatoires présentant un comporte¬ 
ment fractal, pour lesquelles les même motifs se retrouvent à 
différentes échelles, comme c’est le cas pour les nuages ou les 
reliefs montagneux. L’outil mathématique derrière cette pro¬ 
priété d’auto similarité est le mouvement brownien fractionnaire, 
modèle stochastique bien connu introduit dans d. Néanmoins, 
l’utilisation des champs browniens fractionnaires, qui sont iso¬ 
tropes par définition, n’est pas adaptée à des données présentant 
de l’anisotropie. D’autres modèles ont été introduits dans la lit¬ 
térature pour y remédier, citons notamment le drap brownien 
fractionnaire O O et le champ brownien fractionnaire aniso¬ 
trope O, satisfaisant une propriété d’anisotropie globale, ou 
dans un autre registre les textures localement parallèles la.La 
conception mathématique de modèles de textures anisotropes 
et leur synthèse, sont des éléments essentiels à l’estimation des 
caractéristiques anisotropes des images. Nous avons introduit 
dans IHl un nouveau modèle. Nous proposons dans cet article 
deux simulations par champs tangents de ce modèle, dont un 
est extrait de |[6l. La brique de base de ce travail s’appuie sur le 
champ brownien fractionnaire dont nous rappelons la définition 


et les propriétés. 

Le champ brownien fractionnaire (CBF) d’indice de Hurst 
H G (0,1), désigné par , est l’unique champ gaussien à 
valeurs réelles vérifiant : 

- presque sûrement 5^(0) = 0, 

- incréments stationnaires : pour tout z G , 

B^{- + z) - B^(z) = B^(-) - B"(0), 

- autosimilarité d’ordre iJ : VA G K*, B^{-), 

- isotropie : pour toute rotation R de M?, o R = B^ ^ 

jT, 

OÙ = désigne l’égalité en loi. 

D’après |[7)|, le CBF peut se définir à partir de sa représentation 
harmonisable : 

où dW est une mesure brownienne complexe et x • ^ désigne le 
produit scalaire sur L’indice de Hurst H est un paramètre 
fondamental du CBF, comme indicateur de rugosité de la tex¬ 
ture. Plus H est grand, plus la texture résultante est lisse. 

Dans le but d’introduire de l’anisotropie dans le précédent 
modèle, Bonami et Estrade El ont remplacé l’indice de Hurst 
H dans O par une fonction dépendant de la direction ^ et ont 
dérivé une nouvelle classe de champs browniens fractionnaires 
anisotropes (CBFA) : 

r _ 1 ^ 

X(x)= / , 1 ^ 11 ,, (2) 





Figure 1 - Champ élémentaire de paramètres H = 0.2, o^o = 
g et a = 10“^ simulé par (a) la méthode de Cholesky, (b) la 
méthode des bandes tournantes. 


Une classe de modèles anisotropes plus large peut être définie 
par 

X(x)= / (e-'«-l)/V2(^)dïy(4), (3) 

OÙ la densité spectrale / est de la forme 

/i/2(0=c(arg^)||^||-'*(“g«)-i. (4) 

avec c et h deux fonctions 7r-périodiques définies sur l’inter¬ 
valle {—7vl2^ TT/2] dont les images satisfont c((— 7r/2, 7t/2]) C 
M+ et/i((—7r/2,7r/2]) C (0,1). Quand c etsont toutes deux 
constantes, on retrouve le CBF d’ordre H = h. 

Pour définir des modèles anisotropes stationnaires possé¬ 
dant une orientation globale on peut fixer h = H dans dU) 
et considérer : 

Cao,«(arg(^)) = l[-a,a](arg(^) - ao), (5) 

pour un certain 0 < a ^ 7t/2. Nous retrouvons ainsi les 
champs élémentaires de lH, qui sont un cas particulier de CBFA. 
Quand a = 7r/2, ce modèle correspond au CBF isotrope usuel 
d’indice de Hurst H, mais dès que 0 < a < 7r/2, le champ 
n’est plus isotrope dans la mesure où les fréquences non nulles 
sont confinées dans le cône [—a + o^o, + <^o]- Plus la largeur 
du cône a décroît vers 0, plus les fréquences sont concentrées 
suivant l’axe dirigé par o^o, et plus la texture obtenue apparaît 
orientée dans la direction orthogonale à ao, conséquence des 
propriétés de la transformée de Fourier. 

Un exemple de champ d’orientation globale o^o = f est donné 
figure [T] (a) dont la réalisation est obtenue à partir de la fonc¬ 
tion de covariance, qui est couramment employée pour simuler 
des champs gaussiens 0. Cependant la complexité élevée de 
ces algorithmes est un réel problème pour produire de grandes 
images de textures. Une méthode efficace a été proposée dans 
m, appelée méthode des bandes tournantes, et est utilisée ici 
pour simuler la texture de la figure [T] (b). 


2 Une nouvelle classe de champs gaus¬ 
siens à orientation prescrite 

On introduit maintenant un nouveau modèle gaussien comme 
une adaptation locale du champ élémentaire défini dans ([S])-®- 


(0) avec h = H. Plus précisément, on définit le champ brow¬ 
nien fractionnaire localement anisotrope (CBFLA) par : 

X(x)= [ (e“-«-l)/i/2(x,Ot^W(^), (6) 

JR2 

avec 

= Cao,a(x,arg4)||4||“'^“\ (7) 

Cao,a(x,arg^) = 1 (arg(^) - a;o(x)), (8) 

ao étant maintenant une fonction différentiable sur Le CB¬ 
FLA est dérivé du CBFA de la même façon que le mouvement 
brownien multifractionnaire (MBM) est dérivé du mouvement 
brownien fractionnaire (MBF), puisque le paramètre d’orienta¬ 
tion ao est remplacé par une fonction dépendant de la position 
X, tout comme l’indice de Hurst H du MBM peut varier sui¬ 
vant la position lITÜl . 

Champ tangent en chaque point. Afin de décrire et simuler 
cette nouvelle classe de champs gaussiens, nous avons besoin 
de la notion de champ tangent, que l’on rappelle brièvement : 
un champ aléatoire X est dit localement asymptotiquement au¬ 
tosimilaire d’ordre H G (0,1) en xq G si pour tout h G 
le champ aléatoire p~^{X (xq+ ph) —X (xq) ) admet en loi une 
limite Fxo quand p 0 (c.f imiia pour une définition plus 
générale). Le champ Fxo ^st alors appelé le champ tangent de 
X en xo. Ainsi, le champ aléatoire X se comporte localement 
comme Fxo quand x ^ xq. Cette notion a été introduite pour 
la première fois dans lITTIl pour décrire le comportement local 
du mouvement brownien multifractionnaire (qui se comporte 
localement comme un MBF). 

Nous avons prouvé dans {13^ que le CBFLA X de © admet 
pour champ tangent : 

yxo(x) = [ - 1)/'/2(xo,4) dW{^) . (9) 

On remarque que ce champ tangent Fxo ^st un champ élémen¬ 
taire Q-dl-®. La simulation d’un CBFLA requiert mainte¬ 
nant de savoir simuler le champ tangent en tout point d’une 
grille discrète fl [0,1]^ avec r = 2^ — 1, /c G N*. 

Deux méthodes sont proposées : 

(i) Simulation des champs tangents par la méthode de 
Cholesky. La covariance d’un champ élémentaire Fxo ^ fiui est à 
accroissements stationnaires, est reliée à son variogramme 
parCov(rxo(x), Vxo(y)) = (x)( y)( x-y), 

avec déterminé par une formule explicite (H . La matrice 
de covariance H associée est de taille et sa décom¬ 

position de Cholesky H = LL^ est de complexité O(r^), ce 
qui est extrêmement coûteux. Une réalisation du champ élé¬ 
mentaire peut être obtenue par Fxo ^ LZ (figure [T] (a)), où 
Z ^ (A/'(0, l))’^ . Une méthode plus fine est proposée dans O. 

(ii) Simulation des champs tangents par bandes tournantes. 

La méthodologie de IH a été adaptée dans ^ . 

- Formulation discrète du champ tangent 
On peut dériver une expression intégrale du variogramme de 
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Figure 2 - (a) Texture obtenue à partir du champ d’orientations V^l. en jaune et H = 0.2, (b) zoom autour du point rouge 
xo = {x, y) montrant localement un champ élémentaire orienté, (c) diagramme illustrant chaque paramètre du modèle CBFLA. 


1^0 par un changement de variable en polaire : 


^'yxo (x) 



Cao,a(xo,0) |x ■ 


u( 6 >)| 2 ^d 6 > ’ 
( 10 ) 


Oùu(0) = (cos0,sin0) etj(H) = 2 Hr( 2 H"sin(H.) - 
L’intégrale (fTÜb est de la forme ' u(û))dû avec 

V 0 = j(II)cao,a(^o, (*)| ' ■ En ignorant le facteur 

7 (iJ)cQ,o^Q,(xo, 0), on reconnaît que vq est le variogrammed’un 
CBF d’ordre H. Par conséquent, Fxo P^^t être vu comme une 
somme de CBF tournant autour de l’origine. En discrétisant 0 
en un ensemble ordonné {Oi)i^i^n de n bandes d’orientation, 
et en notant {Xi)i^i^n les largeurs de bandes associées = 
Oi-^i — Oi, les champs à bandes tournantes s’écrivent 


E'xo(x) ^ Y^AiC„o,a(xo,0i)Bf (x-u(0i)), (11) 

2=1 


où les sont n CBF d’ordre H. Cette version discrète est 
une bonne approximation dès lors que max Xi ^ e petit. 

i 

- Simulation le long de bandes particulières 

Suivant il, on choisit Si tel que tan( 6 >i) = avec Pi^qi G Z, 

1 Bf f — cos Oi ^ — sin) ; 0 ^ ki,k 2 ^ r 
[\r r J 

COs 6 >A fr^Hn T ^ 

- {B^{kiqi + k 2 Pi); 0 ^ ki,k 2 ^ r} 

rqi J 

et B^(yi • u( 6 >i)) peut ainsi être généré en utilisant l’algorithme 
rapide de Perrin et al ma sur une grille régulière. 

- Choix dynamique des bandes discrètes 

En pratique, le choix des bandes d’orientation {Oi)i^i^n est 
gouverné par le coût de calcul global des B![^ , au travers de la 
programmation dynamique (H . 

Simulation du CBFLA. Comme observé dans ca pour 
le MBM, un champ gaussien peut être simulé à partir de ses 
champs tangents. Le CBELA se comportant localement comme 



Algorithme 1 Simulation du CBELA 

Entrées: r = 2^ — 1, H, ao, a, e 

Sortie: X CBELA de taille (r + 1) x (r + 1) 

F {Pi^Qi)i^i^n ^ ChoixBandesDynamique(r,e) 

2: Calculer et trier les angles {Oi)i^i^ri • Oi ^ atan2(pi, qi) 
3: Calculer les largeurs de bandes • A^ ^ Oi^i—Oi 

4: Générer n CBE : Bf ^ circMBE(r(|pi| + \qi\),H) 

5: Initialisation : X 0 

6: Pour tout {kl, k2) faire 
7: Pour i = 1 à n faire 

8 : ÜOi ^ ^XMH)Ca„a{ikl,k2),0i) 

9 : X{ki, k2) ^ X{ki, k2) + (^iBf {kiqi + k2Pi) 

10: Fin Pour 

11: Fin Pour 


son champ tangent, en chaque pixel xq, on affecte X(xo) = 
^xo (xq) obtenu par la méthode de Cholesky (i) ou des bandes 
tournantes (ii). Le pseudo-code est donné dans l’algorithme [B 
L’étape de pré-traitement (instructions 1,2,3,4 ), laquelle ne dé¬ 
pend pas des orientations locales prescrites, est exécutée une 
fois pour toutes. Le reste de l’algorithme a une complexité li¬ 
néaire en 0(r^ logn). En effet, en chaque point {ki, /C 2 ), une 
bande tournante Oi intervient dans le calcul de X{ki,k 2 ) si et 
seulement si elle appartient au cône CaQ^oi{{ki-,k 2 )^ 0 i) 7 ^ 0 , 
soit \ 0 i — aii{{ki,k 2 ))\ ^ a. Ce tableau Oi étant trié, un tel 
indice i peut être identifié par une recherche dichotomique. 


3 Synthèse de textures orientées 

Les paramètres utilisés dans les simulations sont r = 255, 
a = 10“^, eie = 10“^. Afin d’éviter les artefacts numériques 
dans (Eli on considère une version régularisée de la fonc¬ 

tion indicatrice c^Q^a^ typiquement une gaussienne. Pour 
constant, on retrouve les résultats de (H (figure [T] (b)), qui est 
conforme à ce que l’on obtient par la méthode de Cholesky (fi- 
gure[T](a)) considérée comme exacte. Toutefois l’emploi d’une 








































Figure 3 - Texture obtenue pour le champ et = 0.2. 
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(a) (b) 

Figure 4 - Texture obtenue pour le champ et = 0.7 

(a) par bandes tournantes : bandes d’artefacts (b) par Cholesky 

fenêtre gaussienne lisse la texture et augmente la régula¬ 
rité. Nous présentons maintenant des réalisations de textures à 
orientation prescrite en chaque point xq, donnée par le champ 
de vecteurs Vxq = u(ao(xo)). La figure [2] (a) est la texture 
obtenue à partir du champ de vecteurs = (cos(—7r/2 + 

^),sin(—7r/2 + ^)). Un zoom autour du point xq (figure [2] 

(b) ) montre que localement un CBFLA se comporte comme 
un champ élémentaire. La figure [2] (c) schématise la fonction 
de densité locale au xq et les différents paramètres. Le champ 

y) présente des bandes d’artefacts plus marquées pour des 
indices de Hurst supérieurs {B. = 0.7), comme le montre la fi¬ 
gure [4] (a) tandis que la simulation par la méthode de Cholesky 
en est exempte figure[4](b). On considère enfin le champ de vec¬ 
teurs, = (cos(sin(2x—l)), sin(sin(2x—l)), dont la tex¬ 

ture correspondante est représentée figure[3] D’autres exemples 
ont été donnés dans m. Comme attendu, les textures obtenues 
avec notre approche présentent un comportement local aniso¬ 
trope, orienté orthogonalement au champ de vecteurs. De plus, 
la simulation d’une texture de taille 256 x 256 prend seulement 
quelques secondes avec la méthode des bandes tournantes, et 
plusieurs heures avec la méthode de Cholesky. 


4 Conclusion 

Nous avons introduit un nouveau modèle stochastique pour 
la simulation de textures à orientation locale prescrite, vu comme 


une version locale du champ brownien fractionnaire anisotrope. 
Nous avons tiré parti de la formulation par champ tangent cou¬ 
plée à la méthode des bandes tournantes, pour proposer un al¬ 
gorithme efficace de simulation de ces textures. En couplant 
avec la méthode de Cholesky on supprime les artefacts mais on 
augmente fortement le temps de calcul. Une extension de notre 
modèle à des indices de Hurst variant localement est en cours 
d’étude ça, et servira de test pour des estimateurs de régula¬ 
rité locale ma et d’orientation locale caEi. 
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